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Abstract. We use a dissipative Gross-Petaevskii equation derived from the Bose- 
Hubbard Hamiltonian to study the effect of the thermal component on the stability 
of a current-carrying superfluid state of a Bose gas in an optical lattice potential. We 
explicitly show that the superfluid state becomes unstable at certain quasi-momcntum 
of the condensate due to a thermal component, which is locked by an optical lattice 
potential. It is shown that this instability coincides with the Landau instability derived 
from the GP equation. 
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1. Introduction 

Recently neutral atomic gases in an optical lattice potential have attracted much interest 
in the condensed matter people. This is because an optical lattice potential is the 
ideal realization of the solid state physics such as the formation of energy bands, Bloch 
oscillations, and Josephson effects, and parameters of the system can be easily tuned 
experimentally. More exciting feature appears in the combination of a Bose condensate 
and optical lattice potentials. These systems show the macroscopic quantum phenomena 
due to the phase coherence of the Bose- Einstein condensate such as the Mott insulator- 
superfluid transition realized in recent experiments jH E] . 

In particular, properties of collective modes of a Bose condensed gas in an 
optical lattice are extensively studied both experimentally [3 IU EJ EH El IE] and 
theoretically (HI EH HUES ESI EH EE] in the last few years. Much of attention is attracted 
to the phenomena of instability of a superfluid state. Within the Gross-Pitaevskii 
(GP) theory in a periodic potential, two different mechanisms of the instability are 
known (HI ^3 El El- One is Landau or energetic instability which occurs when 
the excitation energy becomes negative. This kind of the instability is general in 
the superfluidity, not restricted in the lattice system. The other is dynamical or 
modulational instability which occurs when the excitation energy becomes imaginary. 
Such an instability typically exists in non-linear equations. It is important to note 
that although the dynamical instability occurs at zero temperature and can be studied 
using the GP equation, the Landau instability generally occurs at finite temperature. 
Therefore one cannot study the microscopic process of the Landau instability by using 
the GP theory. 

In experiments at finite temperatures such as Refs. jUEHHJ, there exists a thermal 
component around the condensate. As mentioned above, these situations cannot be 
investigated by the usual GP equation since it does not include a thermal component. 
In recent literature about the theory of a Bose gas in an optical lattice, however, detailed 
work considering a role of non-condensed atoms in the instability of the superfluidity 
is lacking. Thus, it is important to include the effect of a thermal component into the 
theory to go beyond the GP theory and give a microscopic physical picture of the Landau 
instability. In the present paper, we use a dissipative GP equation, which includes the 
effect of a thermal component, and explicitly show that non-condensed atoms locked by 
an optical lattice potential induce an instability of a current-carrying superfluid state. 
Our argument clarifies the physical picture of the Landau instability of a Bose gas in an 
optical lattice, the physical process of which is usually ambiguous in general discussions. 

The present paper is organized as follows. In Sec. El we give a dissipative GP 
equation derived from the Bose-Hubbard Hamiltonian in an optical lattice potential, 
which includes a new term expressing collisions between a condensed and non-condensed 
atom. From the dissipative GP equation, we derive damped Bogoliubov equations in 
which an imaginary term due to a thermal component exists and then solve these 
equations to obtain the excitation energies. This procedure is done for a current- 
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free state in Sec. 13.11 and a current- carrying state in Sec. 13.21 We show that due 
to collisions between condensed and non-condensed atoms, the damping (Sec l3.1|) and 
instability (Sec 13.2)1 appear for a current-free and current-carrying state, respectively. In 
Sec. 13.21 we also discuss the instability obtained in the present paper and the relation to 
the previous works about the instability of the usual GP theory. In Sec. EJ we summarize 
our results. 



2. A dissipative Gross-Pitaevskii equation in an optical lattice 

In this section, we derive an equation of motion for a Bose condensate order parameter 
which include the effect of a noncondensed atoms. 

As a model of Bose atoms in a one-dimensional optical lattice, we use the following 
Bose-Hubbard Hamiltonian : 

w = -JY, 0&+i + b l^) + ^E 6 M ft A, (i) 

i i 

where hi and b\ are annihilation and creation operators for Bose atoms at lattice site i, 
respectively. J and U are hopping parameter and onsite interaction, respectively. From 
the Bose-Hubbard Hamiltonian (JTJ), one can derive the dissipative GP equation jTHUHl 

rum] 

ih^Mt) = -J + + {U [<(*) + 2fH{t)] - ihR(Tit)} (2) 

where = is an order parameter of a condensate. In Eq.(J2J, we defined the 
dissipative term R as 

^) = l{^r-) 2 E Hkc-h-h + h) 



ft V2iV S i te / fcljfc2jfe3 



xSi e1{t) - Eiih, t) - Etfa, t) + Ei(k 3 , t) 

x {i + fiih^yii + f^t^f^t) 



(3) 



-/i(M)/i(*2,*){l + /i(A;3,*)} 

The condensate momentum k c and the condensate energy ef(t) will be defined later. 
The noncondensate atoms are treated in the Hartree-Fock (HF) approximation : 

Ei(kj) = -e(k J ) + 2U{ n c(t)+n l (t)}, (j = 1,2,3), (4) 

where e(kj) = 2 J cos(kja) with a being a lattice constant. The dissipative term 
given by Eq. (jHJ) describes collisions between the condensate and the noncondensate 
atoms rjTj. Here, we used the following two approximations in deriving the dissipative 
GP equation (J2J). First, we use the Hartree-Fock approximation for the excitation 
energies of the noncondensed atoms. This Hartree-Fock approximation used in Eq.(TJ]) 
is valid only when Un c /J 1. In this regime, there is no difference between the 
Bogoliubov and Hartree-Fock approximation [21] . Second, we assume the non-condensed 
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atoms occupies only the first band. This assumption is valid when the band-gap between 
the first and second band is large compared to the thermal energy, i.e., 2/c#T <C sEr, 
where s is the strength of the optical lattice potential and Er is the recoil energy. 

The dissipative term Ri(t) is related to the collision term C\ 2 in the Boltzmann 
equation for the noncondensate atoms which describes the dynamics of a 

nonequilibrium distribution function. In the present study, we consider the situation 
where the thermal cloud is locked by an optical lattice potential and thus we assume that 
the thermal cloud is sufficiently close to static equilibrium. Instead of a nonequilibrium 
distribution function, we use the Bose distribution function 

fi(k,t)~f(k) 

= exp {jl [-e(fc) + 2U (n c0 + h°) - fa]} - 1 ' ® 

where n c0 and n° are the equilibrium densities of the condensate and the noncondensate, 
respectively and Jxq is a chemical potential of the non-condensed atoms in equilibrium. 
Denoting the dissipative term in which / is displaced by f° in Eq. (J3J) as i?°(r^) , the 
dissipative GP equation (J2J) becomes 



fo^Mt) = -J[*i+i(t) 



{U [<(*) + 2n,] - ihR°(t)} 



(6) 



Using the equilibrium Bose distribution one finds that the dissipative term is 
riven by [221 123 H3 

-j9[A-«?W] 



n(t) 



(7) 



where we defined the collision time 



1 / U 



nit) h \2K 



site 



) 2 E Hkc 

/ 7,- 7„_ J„„ 



k\ - k 2 + k 3 



x5 (et(t) - Eiih, t) - E t {k 2 , t) + Ei{k 3 , t)) 



xf(h)f°(k 2 ) 



Writing the order parameter ^i(t) in terms of the amplitude and phase, 



(8) 



(9) 



we rewrite the dissipative GP equation (J2J) as 
dt ' 



2JJnt{t)n c i+1 {t) sin (6 i+1 (t) - 4 (t)) 



+2J^jn c i {t)n c i _ 1 {t)sm (6 t (t) - 0i_i(t)) + hn^R , 



(10) 



d0j{t) 
1 8t 



J 



n 



i+1 



(t) 



cos (8 i+1 (t) - 8i(t)) 



cos (8i(t) -6i-i(t)) 
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The right hand side of Eq. gives the energy of the condensed gas, i.e., 

m 



j 



cos (9 i+1 (t) - 9i(t)) + 



nU(t) 



cos (9i(t) - 9i-x(t)) 



\ <{t) 

' -[/«(*) + 2n°). (12) 

In static thermal equilibrium, this energy of the condensate atom reduces to the chemical 
potential of the condensate n° c . 



3. Damped Bogoliubov equation for a current-free and current-carrying 
state 

In this section, we derive Bogoliubov equations from the dissipative GP equation (J2J) 
following Williams and Griffin [221 • First, we derive Bogoliubov equations for a current- 
free state which include a damping term due to a thermal component, which we call 
damped Bogoliubov equations. Second, we generalize the damped Bogoliubov equations 
for a current-carrying state. 



3.1. Damped Bogoliubov equation for a current- free state and damping of the excitation 

In this section, we concentrate on the Bose condensate atoms at rest, i.e., k c = 0. If 
the gas is weakly disturbed from equilibrium, one can consider linearized solution by 
writing 

* i {t) = e- i M n [V + 6V i (t)], (13) 

where /i c0 is a chemical potential of a condensate in static equilibrium. Substituting 
Eq. fjlHj) into the dissipative GP equation (jHJ) and performing some algebra, one obtains 
the linearized equation of motion for a condensate fluctuation <5\l/(£) : 

ih^5Vi(t) = -J [5V i+1 (t) + ^i-i(t)] [2U( 



cO ,,-Ch 



cO 



6Vi{t) + Un c0 5^*(t) 



.Tif3n 



cO 



J 



6V i+1 (t) + 6V* +1 (t) + S%^(t) + S^U(t) 



2n c0 
-2 + 



(14) 



where we defined 



1 / U 



t° h \2N <t 



J2 SiEikJ+EikJ-Eih)) 

x5(k 2 + h- h) fikjfih) {l + f°(h)} . 
In order to solve Eq.([14|). we expand as 



u ^ e i(<iia/ft-u q t) _|_ v * Q-^qla/K-Uqt) 



(15) 
(16) 
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and substitute Eq. (|TBj) into Eq.flTJJ). We then obtain the damped Bogoliubov equations 

h/3 



huJ q U q 



(17) 



u 



<[■ 



where 



A(q) = - 2Jcos(qa) + 2U n cU + n u - fj! 



,c0 



4 J sin' 



qa 



Uri 



cO 



(19) 



Solving the damped Bogoliubov equations (fP7|) and (JTHJ) to first order in l/r°, we find 
that the excitation energy is given by 

2~ 



fru;„ = 2 sin ( ^ ) \ l2UJn c0 + 4 J 2 sin 2 f ^ 



cO 



(20) 



One can see that the first term in this expression (|20|) shows a usual Bogoliubov 
spectrum. For small q, the spectrum ()20j) reveals a phonon-like collective mode. On the 
other hand, the second term in Eq. (j2(J|) works as the damping to the collective mode 
in a condensate since this term is always positive. This type of damping arises from 
the collisions between condensed and non-condensed atoms [221 12H| and exists even if 
non-condensed atoms are in static equilibrium. We note that there exists the Landau 
damping which originates from the interaction between a condensate collective mode 
and the excitations of the thermal component [2J- To include the Landau damping, 
one has to consider the dynamics of noncondensed atoms. 



3.2. Damped Bogoliubov equation for a current- carrying state and instability of the 
superfluidity 

In this section, we derive damped Bogoliubov equations for a current-carrying state, 
i.e., a condensate has a quasi-momentum k c ^ and discuss the instability of the 
superfluidity. 

Instead of Eq. (jl3|) . we substitute the following ansatz 

= e ikJa/h-„°°t/h) ^ q + gq^y ^ ( 21 ) 

into Eq.(J2J). Then one obtains an equation of motion for the fluctuation of the order 
parameter : 
d 

ih—5%(t) = -J 



+ 



e ik ' a / n 5V i+1 (t) + e- ikca/h 5^/ l . 1 (t)] + Un c0 S^*(t) 
2U{n c0 + n°) - fif] 6Vi(t) - ih^ Q 5R°{t), 

cO i n^,0\ 



(22) 



where pj£ = — 2J cos(k c a) + U(n + 2n ) and 

8R°(t) = -4^diff(*) 

g f J 
r° I 2v / ^o 
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cos(k c a) 

5V l+1 (t) - 5V* l+1 (t) - + 8*1$)] sm(k c a) 



-2 (**,(*) + *¥?(*)) 
J 



2iVn c0 

-Uy/n^[S%(t) + S^(t)]^. 
Similarly to the current-free case, we substitute 



5%(t) = E 



i(qla/h-w kq t 



-i(qla/h—uifct\ 



into Eq. (j22j) . We then obtain the equations similar to Eqs. (jl7|) and 

M 7, 



flUJkqUq 



A(k c + q)- t^A(k c - q) 



Un 



i—^A{k c + q) 



%8 ~ 

A(k c -q)+i^A(k c + q) 



v„ + 



Un c0 + i^A(k c -q) 



W 7, 



2r° 



u 



(23) 
(24) 

(25) 
(26) 



(27) 



- hu kq V q 

where 

A(k c ± q) = Un c0 — 2 J cos {(k c ± q)a) — cos(k c a) 

Using the damped Bogoliubov equations for the current- carrying state ()25|) and 
(|26|) . one obtains excitation spectra. First, we derive the excitation spectra from Eq.([25|) 
and ()26|l when the imaginary terms are neglected. After some algebra, one obtains 

huJkq = 2Jsin(fc c a) sin(ga) 



±2 sin 



qa 

T 



' 2UJn c0 cos(k c a) + 4 J 2 cos 2 (A; c a) sin 2 



qa 

T 



(28) 



From this expression of the excitation spectra, one can see that the instability occurs 
when the excitation energy becomes negative, or when it develops the imaginary part, 
i.e., k c > 7r/2. One can identify the origin of the first instability as the Landau 
instability, and the second one as the dynamical instability, which were discussed in 
detail in Ref. PEED]. 

Second, we calculate the excitation spectra including the imaginary terms. We 
obtain the following expressions : 

huJkq = 2Jsin(fc c a) sin(ga) 



±2 sin 
.2nd 



qa 

T 



'2UJn c0 cos(k c a) + 4J 2 cos 2 (fc c a) sin 2 



qa 

T 



2 J cos(k c a) sin 2 (%r) + J 



, OL 



a + 1 



sin(/c c a) sin(ga) + Un 



cO 



(29) 



where 



4J 



a a 
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Stability Phase Diagram : Urf/J = D.1 Stability Phase Diagram : lilt"/ J = 015 




a 112 oa tiA ub i n H p.t 



kclq s : Condensate Quasi-Momentum kc/qg : Condensate Quasi-Momentum 

Figure 1. Stability phase diagram. The horizontal axis is the quasi-momentum of 
the condensate and the vertical axis is the quasi-momentum of the excitation, qg is 
the Bragg wave vector defined by qs = ft/a. Colored regions are unstable : Blue area 
is the dynamically unstable region and both red and blue areas are the Landau and 
dissipative unstable region. 

Eq. (l29J) is the main result of the present paper. From Eq. (j29|) . one can find a new 
origin of the instability, which we call "dissipative instability". When the sign of 
the third term of Eq. (j29|) is changed at a certain quasi-momentum k c , the excitation 
energy increases exponentially as time increases, which means that the superfluid state 
becomes unstable. This instability directly comes from the collisional term between 
the condensed and non-condensed atoms, and thus arises due to the mutual interaction 
between the condensed and noncondensed atoms. The stable phase diagrams in k c -q 
plane are shown in Fig^ From Fig^ one can see that the critical line of the Landau 
instability, which is derived from the GP equation coincides with the critical line 

of the "dissipative instability" . This is physically reasonable because in experiments the 
thermal component should be playing a role of the obstacle that induces the Landau 
instability. Here, we have explicitly shown that the collisions between the condensed 
and non-condensed atoms locked by an optical lattice potential indeed make the system 
unstable in the sense as the Landau instability. We note that this kind of instability 
is analogous to the instability of a condensate induced by a rotating thermal cloud 
discussed by Williams et al. [24. 

4. Summary 

In the present paper, we have studied the instability of the current-carrying superfluid 
state at finite temperature in an optical lattice potential. Using the dissipative GP 
equation derived from the Bose-Hubbard Hamiltonian, we have obtained the damped 
Bogoliubov equations and calculated the excitation spectrum including the imaginary 
part due to a thermal component. By analyzing the excitation spectrum, we have 
explicitly shown, for the first time, that the mutual friction between the condensate and 
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the thermal component locked by an optical lattice potential leads to the instability at a 
critical condensate momentum, which is equivalent to the Landau instability criterion. 
This work gives the clear physical picture about the Landau instability of a Bose gas in 
an optical lattice at finite temperature in the actual experimental situation. 

In the present work, we have only investigated the stable phase diagram focusing 
on the effect of the thermal component. It will be important and interesting to compare 
our analysis with the experimental data performed recently [H|. In this connection, it 
will be interesting to investigate the temperature dependence of the imaginary term in 
the excitation spectrum and calculate the growing time of the excitations. 
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